Available  online  at  www.sciencedirect.com 


ELSEVIER 


%V  ScienceDirect 


Journal  of  Power  Sources  161  (2006)  1041-1055 


JOURNAL  OF 


www.elsevier.com /locate /jpowsour 


Optimization  of  a  PEM  fuel  cell  system  based  on  empirical  data 
and  a  generalized  electrochemical  semi-empirical  model 

J.  Wishart,  Z.  Dong*,  M.  Secanell 

Department  of  Mechanical  Engineering  and  Institute  for  Integrated  Energy  Systems,  University  of  Victoria,  Victoria,  BC,  Canada  V8W  3  P6 

Received  10  March  2006;  received  in  revised  form  16  May  2006;  accepted  16  May  2006 

Available  online  17  August  2006 


Abstract 

A  systematic  method  to  obtain  the  optimal  operating  conditions  of  a  fuel  cell  system  is  presented.  This  method  is  based  on  the  coupling  of  a 
semi-empirical  fuel  cell  stack  model  and  an  associated  balance  of  plant  (BOP)  model  with  an  optimization  algorithm  in  order  to  efficiently  explore 
the  range  of  possible  operating  conditions.  The  approach  described  in  the  paper  to  obtain  optimal  operating  conditions  is  applied  to  a  fuel  cell 
system  designed  to  operate  in  two  different  applications:  automotive  and  stationary.  In  both  cases,  the  application  of  this  methodology  results  in  a 
set  of  optimal  operating  conditions  that  yields  large  improvements  in  the  system  performance.  The  optimization  problem  is  solved  for  two  different 
performance  objectives:  maximization  of  net  system  power  and  maximization  of  system  exergetic  efficiency. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

As  a  promising  technology  that  may  successfully  supersede 
the  combustion  of  fossil  fuels  as  the  dominant  method  of  energy 
conversion,  hydrogen  fuel  cells  are  studied  worldwide  with  an 
aim  to  improve  the  power  output,  lower  the  cost  and  extend  the 
life  of  operation  for  widespread  applications. 

Among  various  types  of  fuel  cells,  the  proton  exchange  mem¬ 
brane  fuel  cell  (PEMFC)  is  arguably  the  fastest-growing  type 
and  the  fuel  cell  that  is  most  likely  to  be  widely  used  in  the  near 
future.  The  modeling  and  optimization  of  a  PEMFC  system,  car¬ 
ried  out  in  this  work,  is  aimed  at  achieving  better  performance 
of  a  given  fuel  cell  system  design. 

1.1.  Stack  modeling 

Modeling  of  real-world  applications  has  been  seen  as  a  useful 
tool  for  decades.  Fuel  cell  and  fuel  cell  system  modeling  is  in  its 
relative  infancy,  but  already  a  significant  amount  of  effort  has 
been  put  forth  to  understand  the  parameters  and  issues  affecting 
the  performance  of  the  fuel  cell  system. 
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There  are  several  types  of  modeling  approaches,  and  the  line 
separating  these  approaches  is  often  blurry.  Essentially  however, 
the  approaches  may  be  classified  as:  (a)  theoretical,  sometimes 
known  as  ‘mechanistic’;  (b)  computational  fuel  cell  dynamics 
(CFCD)  simulation;  (c)  semi-empirical;  (d)  empirical  depend¬ 
ing  on  the  level  of  modeling  sophistication.  Each  approach  has 
advantages  and  disadvantages,  as  discussed  in  [1]. 

In  this  work,  since  the  main  interest  is  in  the  development 
of  a  methodology  to  obtain  the  optimal  fuel  cell  system  oper¬ 
ating  conditions  and  not  on  the  fuel  cell  design  itself,  a  simple 
but  computationally  inexpensive  model  is  used.  In  particular, 
the  semi-empirical  approach  based  on  work  by  researchers  at 
the  Royal  Military  College  (RMC)  is  chosen  [2-5].  The  RMC 
model  is  robust  and  flexible  and  is  primarily  steady-state  but 
also  has  simplified  transient  aspects.  The  RMC  model  has  been 
used  by  members  of  industry,  where  Ballard  Power  Systems 
is  the  most  notable  example  [6],  and  by  other  research  groups 
[7]  that  have  incorporated  early  versions  of  the  model  and  used 
experimentation  to  validate  its  veracity.  The  model  has  achieved 
an  accuracy  and  adaptability  that  has  allowed  it  to  be  deemed 
functional  for  the  current  research. 

In  order  to  make  this  semi-empirical  model  as  accurate  as 
possible,  over  the  past  1 0  years,  considerable  work  has  been  done 
to  acquire  the  empirical  model  parameters  of  the  benchmark 
Ballard  MK4  and  MK5  stacks.  Extensive  tests  have  also  been 
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Nomenclature 


Ph2o 


/^Interface 

co2 


exCH 

oxke 

expE 

-An  ax 


^Nernst 

/(x,  i) 

F 

h 

^conv 


^active  active  cell  area  (cm2) 

Acen  radiative  stack  area  (cm2) 

Agnds  area  of  fuel  cell  stack  ends  (cm2) 

Asides  area  of  fuel  cell  stack  sides  (cm2) 

At  total  stack  cross-sectional  area  (cm2) 
cp ^  specific  heat  of  air  at  constant  pressure 

(JK-1  kg-1) 

specific  heat  of  water  at  constant  pressure 
(4128  JK"1  kg"1) 

^interface  concentration  of  oxygen  gas  at  the  surface  of  the 
catalyst  at  the  cathode 
ex  exergy  (kJ  mol  - 1 ) 

excH  chemical  exergy  (kJ  mol-1) 

exKE  kinetic  energy  exergy  (kJ  mol-1) 
expE  potential  energy  exergy  (kJ  mol-1) 

£max  maximum  voltage  obtained  from  converting 

enthalpy  of  formation  of  hydrogen  into  electrical 
energy  (1.48  V) 

^Nernst  Nernst  voltage  (V) 
fix,  i )  objective  function 

F  Faraday  constant  (96,485  C  electron- 1 ) 
h  specific  enthalpy  (kJ  mol-1) 

hcony  heat  transfer  coefficient  of  convection 
(Wm  2  K-4) 

/zConv, bottom  heat  transfer  coefficient  of  convection  of  fuel 
cell  channel  bottom  (W  m-2  K-4) 

^conv,top  heat  transfer  coefficient  of  convection  of  fuel  cell 
channel  top  (W  m-2  K-4) 

/zprodj  enthalpy  of  product  j  (kJ  mol-1) 

^  react,/  enthalpy  of  reactant  i  (kJ  mol-1) 
ho  enthalpy  of  reference  environmental  (restricted) 

state  (kJ  mol-1) 

i  fuel  cell  current  density  (A  cm-2) 

I  fuel  cell  current  (A) 

/l  current  at  which  hydrogen  consumption  is  equal 
to  hydrogen  supply  (A) 

/vap  latent  heat  of  vaporization  (J  kg  - 1 ) 
ma |r  mass  flow  of  air  (kg  s-1) 

fi^cool,H2o  mass  flow  of  water  required  for  cooling  (kg  s  - 1 ) 
rhu2  mass  flow  of  hydrogen  gas  (kg  s-1) 
fi^H2o,heat  exchanger  mass  flow  of  water  in  heat  exchanger 

(kgs  A 

mH2o,in  air  mass  flow  of  water  required  for  air  stream 
humidification  (kg  s-1) 

fi^H2o,in  h2  mass  flow  of  water  required  for  hydrogen 
stream  humidification  (kg  s-1) 
fi^H2o, produced  rate  of  water  produced  at  cathode  (kg  s-1) 
Ma u-  molar  mass  of  air  (28.97  x  10-3  kg  mol-1) 

Mu2  molar  mass  of  hydrogen  (2.016xl0-3  kg  mol-1) 
Mh2o  molar  mass  of  water  (18.018  x  10-3  kg  mol-1) 
nCQ n  number  of  fuel  cells  in  the  stack 

•  I 

Ah2  molar  flow  of  hydrogen  gas  at  inlet  (mol  s  ) 


"prodj 
^  react,/ 

ho 


^interface 

p  o2 


nsat 

^h2o 


P  air, atm 
^cell 


^cell 

Vh2 


^H2,out  molar  flow  of  hydrogen  gas  at  outlet  (mol  s 


•  I 

No2  molar  flow  of  oxygen  gas  at  inlet  (mol  s  ) 

•  t 

Aprod  j  molar  flow  of  product  j  (mol  s  ) 

Areact ,/  molar  flow  of  reactant  i  (mol  s-1) 

^interface  pressure  0f  hydrogen  gas  at  surface  of  the  catalyst 

at  the  anode  (bar) 

^mterface  pressure  0f  oxygen  gas  at  surface  of  the  catalyst 
at  the  cathode  (bar) 

P^2q  saturation  pressure  of  water  vapour  at  a  given 
temperature  (bar) 

^air,atm  atmospheric  pressure  (101.325  kPa) 

Pceii  fuel  cell  operating  pressure  (bar) 

^drop, coolant  pressure  drop  of  coolant  through  cooling  loop 
(bar) 

^drop, humid  pressure  drop  of  water  through  humidification 
stream  (bar) 

em  membrane- specific  resistivity  for  the  flow  of 

hydrated  protons  (Q  cm) 

R  universal  gas  constant  (8.3145  J  mol-1  K-1) 
^electronic  0hmic  electronic  overvoltage  resistance 

(Q  cm2) 

^protonic  ohmic  protonic  overvoltage  resistance  (Q  cm2) 
Gair, total  heat  lost  to  air  stream  and  surroundings  (W) 
Gconv, total  heat  lost  to  air  stream  due  to  convection  (W) 

Q conv,top+bottom  heat  lost  to  convection  at  the  fuel  cell 
channel  top  and  bottom  (W) 

Gconv, walls  heat  lost  to  convection  at  the  fuel  cell  channel 
walls  (W) 

Gevap  heat  lost  to  air  stream  due  to  evaporation  (W) 
Grad  heat  lost  to  surroundings  due  to  radiation  (W) 

G total  heat  produced  due  to  electrochemical  reaction 

(W) 

s  specific  entropy  (kJ  mol-1  K-1) 

so  entropy  of  reference  environmental  (restricted) 

state  (kJ  mol-1  K-1) 
t  thickness  of  membrane  (cm) 

Teen  fuel  cell  operating  temperature  (K) 

7cooi,H2o,in  temperature  of  coolant  water  at  inlet  (K) 

7o  ambient  air  temperature  (K) 

Vceii  output  voltage  of  single  fuel  cell  (V) 

^compressor  power  consumed  by  the  compressor  (W) 
Wcool,pump  power  consumed  by  the  coolant  pump  (W) 
fkhumid,pump  power  consumed  by  the  humidifier  pump  (W) 
fkgross, stack  gross  power  output  of  fuel  cell  stack  (W) 
wj  weighting  of  objective  j 

Wfix,  i )  net  power  at  current  density  ij  for  objective  j  (W) 
Wnet, parasitic  net  parasitic  power  lost  (W) 

Wnet, system  net  power  output  of  fuel  cell  system  (W) 
x  design  variable  vector 

Xj  molar  fraction  of  substance  j 

Xmax  upper  limit  vector  of  design  variable 

xmin  lower  limit  vector  of  design  variable 

xh2o  molar  fraction  of  water  in  a  gas  stream  for  a  given 
temperature 

^otheTgases  m°lar  fraction  of  gases  (apart  from  oxygen)  in 
the  air  stream 


Gevap 

Grad 

G  total 


^cell 


Amax 

Xmin 

^.sat 

xh2o 


.channel 
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^othergases  m°lar  fraction  of  gases  (apart  from  oxygen)  at 
the  inlet 

xotoUgases  m°lar  fraction  of  gases  (apart  from  oxygen)  at 
the  outlet 


Greek  symbols 

Pi 


empirical  term  in  activation  overvoltage  equation 
(V) 

empirical  term  in  activation  overvoltage  equation 
(VK-1) 

empirical  term  in  activation  overvoltage  equation 
(VK"1) 

empirical  term  in  activation  overvoltage  equation 
(VK"1) 

isentropic  compression  ratio 
change  in  enthalpy  (kJ) 
activation  overvoltage  (V) 
concentration  overvoltage  (V) 
efficiency  of  model  pump  motors 
ohmic  overvoltage  (V) 
efficiency  of  model  pumps 
system, energy  system  energy  (1st  law)  efficiency 
tj system, exergetic  system  exergetic  (2nd  law)  efficiency 
^ohmic°mC  °hrmc  electronic  overvoltage 

%hmic  ohmic  protonic  overvoltage 

Aair  air  stoichiometry  (=  oxygen  stoichiometry) 

Ah7  hydrogen  stoichiometry 

^membrane  empirical  parameter  that  describes  membrane 
conditions 

oxygen  stoichiometry  (=  air  stoichiometry) 
chemical  potential  of  substance  j  at  reference 
environmental  (restricted)  state  (kJmol-1) 
chemical  potential  of  substance  j  at  reference  dead 
environmental  (unrestricted)  state  (kJmol-1) 
Ph2o  density  of  water  (998  kg  m3) 


ft 

ft 

ft 

y 

AH 

7?  act 
Oconc 
0  motor 
0  ohmic 
77pump 


^o2 

fto 

M/00 


carried  out  in  the  author’s  group  to  acquire  the  model  parameters 
for  the  new  close-to-ambient  pressure  PEM  fuel  stack  designs, 
including  the  tri-flow,  external-manifold,  radiator  stack  (TERS), 
developed  by  the  group  [8-10],  and  the  PC6-1200  stack  from 
Palcan  Power  Systems  [11,12]. 

A  PC6-1200  stack  has  been  installed  as  part  of  a  power  sys¬ 
tem  in  a  low-speed  fuel  cell  hybrid  electric  vehicle  (LSFCHEV). 
This  LSFCHEV  will  undergo  performance  testing  using  a 
dynamometer  platform  that  will  provide  indirect  system  data, 
and  the  system  itself  will  be  tested  separately  to  complement 
the  data  set  and  measure  the  accuracy  of  the  system  model.  Sub¬ 
sequent  publications  will  include  discussion  and  results  from 
these  tests.  The  data  obtained  from  the  stack  testing  of  the 
three  aforementioned  stacks  have  been  used  to  determine  the 
model  parameters  associated  with  each  particular  stack,  and  the 
results  are  presented  in  Section  3.1.  The  current  work  is  the 
first  step  in  developing  a  virtual  prototype  of  the  LSFCHEV 
that  will  undergo  extensive  simulation  and  sensitivity  analysis 


before  the  physical  prototype  is  tested  to  validate  the  simulation 
results. 

7.2.  System  optimization 

Given  the  complexity  of  fuel  cell  modeling,  fuel  cell  and 
fuel  cell  system  optimization  presents  a  real  technical  challenge. 
Nevertheless,  some  laudable  attempts  at  fuel  cell  optimization 
can  be  found  in  the  literature.  Grujicic  and  Chittajallu  [13]  used 
a  two-dimensional  computational  fuel  cell  dynamics  model  of 
a  single  fuel  cell  to  optimize  the  electric  current  per  fuel  cell 
width  at  a  cell  voltage  of  0.7  V.  In  the  optimization,  sequential 
quadratic  programming  was  used  to  obtain  the  operational  and 
geometric  parameters  for  achieving  the  maximum  electric  cur¬ 
rent,  including  air  inlet  pressures  and  cathode  thickness,  cathode 
length  for  each  shoulder  segment  of  flow  channel,  and  fraction  of 
cathode  length  associated  with  the  flow  channel.  More  recently, 
Mawardi  et  al.  [14]  employed  a  numerical  one-dimensional  fuel 
cell  model  to  optimize  the  power  density  with  respect  to  several 
operating  parameters.  The  optimization  allowed  for  quantitative 
discussion  of  the  effects  of  membrane  and  electrode  thicknesses 
and  CO  concentration  on  the  values  of  the  optimum  operating 
conditions.  At  the  fuel  cell  system  level,  however,  the  only  two 
attempts  at  optimization  were  done  by  Xue  and  Dong  [6],  and 
Wang  et  al.  [15].  Xue  and  Dong  used  a  semi-empirical  model 
of  the  Ballard  Mark  IV  fuel  cell  and  models  for  the  auxiliary 
systems  to  create  a  model  of  the  fuel  cell  system.  Using  this 
model  and  numerical  optimization,  the  optimal  active  stack  area 
and  air  stoichiometric  ratio  was  obtained  to  maximize  net  power 
output,  and,  at  the  same  time,  minimized  production  costs.  Wang 
and  Dong  performed  an  optimization  on  a  fuel  cell  system  as  a 
demonstration  of  a  novel  optimization  algorithm  [8,15].  The 
optimization  results  showed  an  improvement  in  system  effi¬ 
ciency,  net  power,  volumetric  and  gravimetric  power  densities, 
as  well  as  in  a  reduction  of  system  costs. 

In  general,  optimization  of  fuel  cell  systems  is  still  a  chal¬ 
lenge  not  only  because  of  the  inaccuracy  of  the  models  but 
because  the  optimization  is  a  highly  non-linear  problem  where 
the  objective  function  is  calculated  using  a  numerical  model 
of  the  fuel  cell  and  fuel  cell  system.  Non-linear  optimization 
involves  the  search  for  a  minimum  of  a  non-linear  objec¬ 
tive  function  subject  to  non-linear  constraints.  It  is  common 
for  these  optimization  problems  to  have  multiple  optima.  At 
present,  two  different  search  approaches  have  emerged  in  the 
area  of  non-linear  design  optimization:  local  methods  and  global 
methods. 

Local  methods  aim  to  obtain  a  local  minimum,  and  they  can¬ 
not  guarantee  that  the  minimum  obtained  is  the  absolute  mini¬ 
mum  for  a  non-unimodal  objective  function  and/or  a  non-convex 
feasible  region.  These  methods  are  usually  first-order  methods. 
Some  of  the  most  popular  local  methods  for  optimization  include 
the  conjugate  gradient  algorithms  and  the  quasi-Newton  meth¬ 
ods  for  unconstrained  optimization  and  the  sequential  linear  and 
quadratic  programming  methods  for  constrained  optimization. 
Although  local  methods  do  not  aim  for  the  global  optima,  sev¬ 
eral  approaches  can  be  used  to  continue  searching  once  a  local 
minimum  has  been  located  to  obtain  the  global  optima,  such  as 


1044 


J.  Wishart  et  al.  /  Journal  of  Power  Sources  161  (2006)  1041-1055 


the  stochastic  based  approaches  of  random  multi-start  methods 
[16,17]  and  ant  colony  searches  [18]. 

Global  methods  aim  to  obtain  the  global  minimum  of  the 
function.  Mostly  based  on  stochastic  procedures,  these  meth¬ 
ods  do  not  need  any  information  about  the  gradient.  Some  of 
the  most  popular  methods  for  global  optimization  are  genetic 
algorithms  (GA),  simulated  annealing  (SA),  tabu  search,  and 
stochastic  programming. 

In  this  work,  one  local  optimization  algorithm  and  two  global 
optimization  algorithms  are  coupled  to  a  fuel  cell  system  model 
in  order  to  explore  the  design  space  and  to  guarantee  a  solution 
that  is  independent  on  the  method  used.  The  three-optimization 
techniques  are:  sequential  quadratic  programming  (SQP),  sim¬ 
ulated  annealing  (SA)  and  genetic  algorithms  (GA).  The  opti¬ 
mization  algorithms  use  operating  parameters  such  as  the  tem¬ 
perature,  pressure  ratios  and  the  stoichiometries  of  the  reactants 
as  design  variables.  These  optimization  routines  are  used  to 
search  for  the  best  operating  parameters  that  lead  to  peak  fuel 
cell  system  performance.  This  work  has  focused  on  obtaining 
the  optimal  operating  conditions  of  the  fuel  cell  system  and  not 
on  obtaining  optimal  physico-chemical  parameters  of  the  fuel 
cell.  These  are  two  different  goals.  The  objective  of  the  current 
work  is  to  show  a  method  to  select,  for  a  given  fuel  cell,  the  best 
operating  parameters.  Optimizing  the  physico-chemical  param¬ 
eters,  while  also  an  extremely  useful  and  necessary  endeavour, 
has  the  goal  of  improving  the  actual  design  of  the  fuel  cell  and 
stack;  future  work  should  combine  the  optimization  of  operating 
and  physico-chemical  parameters  resulting  in  a  truly  optimized 
fuel  cell  system. 

2.  Fuel  cell  system  model 

The  PEMFC  stack  model  implemented  in  this  work  is  based 
on  the  work  at  the  Royal  Military  College  of  Canada  by  Amphlett 


et  al.  [2,3],  Mann  et  al.  [4]  and  Fowler  et  al.  [5].  The  model, 
known  as  the  generalized  steady- state  electrochemical  model 
(GSSEM),  is  zero-dimensional,  semi-empirical,  isothermal  and 
static  in  nature;  thus,  the  parameters  of  the  equations  are  deter¬ 
mined  experimentally  to  provide  the  time-independent  polari¬ 
sation  curves,  power  curves  and  system  efficiencies  at  various 
operating  conditions.  The  stack  model  is  integrated  into  the 
larger  fuel  cell  system  model  with  its  constituent  BOP  compo¬ 
nents,  and  the  resultant  model  is  able  to  calculate  the  net  output 
power  of  the  system.  A  schematic  representation  of  the  fuel  cell 
system  is  shown  in  Fig.  1. 

2.7.  Fuel  cell  model 

The  fuel  cell  model  is  for  a  PEMFC,  which  uses  the  following 
electrochemical  reaction: 

H2  +  ^02  H20  (1) 

where  the  reaction  is  exothermic. 

The  voltage  of  the  fuel  cell,  Vceii,  is  modeled  as  [4] 

Tcell  —  T^Nernst  T  Oact  T  ^?ohmic  T-  Oconc  (2) 

where  FNernst  is  the  Nernst  equation,  which  is  an  expression  for 
the  electromotive  force  (EMF)  for  given  product  and  reactant 
activities;  rjact  the  activation  overvoltage,  which  is  the  amount 
of  voltage  used  to  drive  the  reaction;  77ohmic  the  ohmic  over¬ 
voltage,  which  is  the  amount  of  voltage  lost  to  the  resistance  to 
electron  flow  in  the  electrodes  and  the  resistance  to  ion  flow  in 
the  electrolyte;  r/conc  is  the  concentration  overvoltage,  which  is 
the  voltage  lost  when  the  concentration  of  reactant  at  the  elec¬ 
trode  is  diminished.  In  Eq.  (2),  the  overpotentials  are  all  positive 
because  they  are  assumed  to  be  losses  of  the  fuel  cell.  This  is  in 
accordance  with  the  derivation  in  Ref.  [2],  where  a  more  detailed 
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Fig.  1.  Fuel  cell  system  schematic  layout. 
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explanation  of  how  these  terms  are  obtained  based  on  electro¬ 
chemical  arguments  is  given.  The  terms  of  Eq.  (2)  are  discussed 
in  the  following  sections  for  the  sake  of  completeness. 

The  output  of  the  fuel  cell  system  will  depend  strongly  on  the 
amount  of  reactants  input  into  the  system;  the  relative  amounts 
that  take  part  in  the  electrochemical  reaction  of  Eq.  (1)  are  dic¬ 
tated  by  the  molar  balance  of  this  equation.  The  molar  mass  flow 
of  H2  is  given  by 


^H2  ft  cell 
IF 


where  A.h2  is  the  stoichiometry  of  the  hydrogen  gas  and  F  is 
Faraday’s  constant.  The  corresponding  expression  for  the  molar 
mass  flow  of  oxygen  is 


/A.  o2  ft  cell 

4 F 


The  oxygen  flow  comes  from  the  ambient  air  brought  into  the 
compressor  such  that  Aair  is  equal  to  Xo2 . 

The  production  of  water  is  both  a  benefit  and  a  hindrance 
to  the  operation  of  the  fuel  cell,  as  a  delicate  balance  must  be 
struck  between  having  enough  moisture  to  maintain  the  humid¬ 
ity  of  the  membrane  while  also  ensuring  that  the  MEA  is  not 
flooded  with  excess  water.  Either  situation  will  result  in  a  reduc¬ 
tion  in  the  efficacy  of  the  fuel  cell  operation.  If  the  membrane 
dries  out,  the  proton  conductivity  and  oxygen  reaction  kinetics 
decrease  temporarily,  but  excessive  dehydration  can  result  in 
permanent  damage  to  the  membrane  and  permanently  impaired 
power  production.  The  model  used  in  this  work  assumes  that 
the  membrane  is  fully  hydrated  without  any  water  management 
issues,  an  assumption  justified  in  other  work  [18].  The  water 
produced  at  the  cathode  is  used  to  externally  humidify  the  reac¬ 
tants,  as  shown  in  Fig.  1,  and  the  equation  for  water  production 
is  given  by  [19] 

Incei\Mii2o 

ft2  H20,  produced  —  0  ^  w) 


where  Mh2o  is  the  molar  mass  of  water  (which  has  a  value  of 
18.018x  10-3  kg  mol-1). 


2.1.1.  Nernst  voltage 

The  Nernst  equation  for  the  reaction  described  above  is  given 
by  [4] 

^Nernst  =  1-229  -  0.85  x  10_3(7’cen  -  298.15)  +  4.3085 

x  10-5rcell(ln(p“e)  +  0.5  ln(/7gf face))  (6) 

where  T^eii  is  the  stack  temperature  (K),  and  pj!j!erface  and 

<rface  are  l^e  hydrogen  and  oxygen  partial  gas  pressures  (bar) 
at  the  surface  of  the  catalyst  at  the  anode  and  cathode,  respec¬ 
tively. 

The  partial  pressures  at  the  catalyst  surface  are  assumed  to 
be  the  same  across  the  entire  cell,  and  are  given  by  [2,6] 


interface 

Pu2 
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1 .653//  Tij}34  channel 
e  e  XH20 


and 

interface 

p  o2 


P cell  [  i 


sat 

h2o 


channel 
•^other  gases 


e0.291//TX] 


where  is  the  molar  fraction  of  water  in  a  gas  stream  at 
saturation  for  a  given  temperature,  v^"rg^ses  the  molar  frac¬ 
tion  of  other  gases  (apart  from  oxygen)  in  the  air  stream,  p^0 
the  saturation  pressure  of  water  vapour  at  a  given  temperature, 
Pceii  the  cell  operating  pressure,  and  i  is  the  current  density 
(A  cm-2).  The  derivation  of  Eqs.  (7)  and  (8)  is  based  on  the 
Maxwell-Stefan  equations,  and  the  interested  reader  is  referred 
to  Ref.  [2]  for  a  more  detailed  explanation.  After  assuming  that 
the  flux  of  water  vapour  and  nitrogen  (in  the  cathode)  are  zero, 
the  Maxwell-Stefan  equation  becomes  a  simple  ordinary  differ¬ 
ential  equation  that  can  be  integrated  analytically.  The  equations 
are  derived  by  solving  this  equation,  substituting  the  value  of 
all  constants  and  approximating  the  diffusion  coefficient  using 
kinetic  theory.  A  detailed  explanation  of  how  the  equations  are 
derived  is  given  in  Ref.  [2]. 

The  molar  fraction  of  water  at  saturation  in  a  gas  stream  for 
a  given  temperature  is  given  by 


_sat 

..sat  _  ftH20 

H2°  "  FF 


where  Pcen  is  the  cell  operating  pressure  and  the  p^o  term  *s 
determined  in  a  fuel  cell  by  the  following  empirical  equation 
[20] 


ln(pfg0)  =  70.434643  - 


7362.6981 

7cell 


+  0.006952085  x  rceU 


-  9.0000  x  In  rceu 


(10) 


The  molar  fraction  of  other  gases  (mostly  nitrogen  gas)  in  the 
air  stream  is  given  by  a  log  mean  average  between  the  molar 
fraction  of  other  gases  in  a  humidified  stream  of  air  at  the  inlet 
and  the  molar  fraction  at  the  outlet  [2,21] 


channel 
•^other  gases 


where 


in, hum  out, hum 

jr  _  Y 

other  gases  other  gases 


ln(v 


in, hum 
other  gases 


out, hum  \ 
other  gases  1 


(ii) 


in, hum  —  —  Tsat  W  0  7Q 

other  gases  'T  •*H20'  V).ly 


and 


out, hum 


i  _  sat 

1  xh2o 


other  gases  ,  +  _  ,  /Xair)(0.2 1  /0.79) 
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03) 


In  Eqs.  (12)  and  (13),  the  0.79  term  refers  to  the  dry  molar 
fraction  of  other  gases  in  air,  while  in  Eq.  (13),  the  0.21  term 
refers  to  the  dry  molar  fraction  of  oxygen  in  air.  The  Aair  term 
denotes  the  stoichiometry  of  the  air  stream. 


2.1.2.  Activation  overvoltage 

The  semi-empirical  equation  for  the  activation  overvoltage  is 
given  by  [4] 

act  =  01  +  02  7ceii  +  03rcen  ln(Cgfrface)  +  04rcell  ln(7) 


(14) 
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Table  1 


Empirical  coefficients  for  PEMFC  stack  model 


Pi 

Pi 

P3 

Pa 

Amphlett  et  al.  [2]  and  Xue  and  Dong  [6]  for  Ballard  Mark  V  stack 

-9.514 

0.00312 

0.0000740 

-0.000187 

Pastula  [20]  and  Wang  and  Dong  [8]  for  an  ambient  air  PEMFC  stack 

-0.944 

0.00354 

0.0000785 

-0.000196 

Mann  et  al.  [4]  and  Fowler  et  al.  [5]  for  Ballard  Mark  IV  and  V  stacks 

-0.948 

a 

0.000076 

-0.000193 

a  ft  =  0.00286  +  0.0002  In  AactiVe  +  (4.3  x  10“5)lnc*v 


where 


^Interface 

co2 


^interface 

P  o2 


5.08  x  106e<-498/7ieii) 


(15) 


and  where  the  expression  for  ^g*erface  is  given  in  Eq.  (8)  and 
the  /3  coefficients  of  Eq.  (14)  are  empirically  determined  for 
each  individual  fuel  cell  stack.  Their  previously  published  values 
are  given  in  Table  1.  The  values  used  during  this  research  are 
contained  in  the  third  entry  of  Table  1 .  The  term  Aactive  in  the 
description  of  is  the  active  area  of  the  fuel  cell  (cm2),  given 
in  Eq.  (18). 


2.1.3.  Ohmic  overvoltage 

The  ohmic  overvoltage  can  be  expressed  in  accordance  with 
Ohm’s  law  as  [4] 


_  electronic  ,  protonic  _  _  •/ oelectronic  ,  r>protonic\ 

'/ohmic  —  '/ohmic  '/ohmic  —  t[n  >  A  ) 

(16) 

where  7?electromc  is  assumed  to  be  a  constant  over  the  operation 
temperature  of  the  PEMFC.  The  electronic  resistance  is  further 
assumed  to  be  inconsequential  in  comparison  to  protonic  resis¬ 
tance,  and  is  thus  ignored.  The  term  ^Protomc  is  known  to  be  a 
complex  function  of  water  content  and  distribution  in  the  mem¬ 
brane,  which  in  turn  is  a  function  of  the  cell  temperature  and 
current.  A  general  expression  for  the  ohmic  resistance  of  the 
electrolyte  is  given  by  [20] 


^protonic  _ 


mt 

Aactive 


(17) 


where  cm  is  the  membrane- specific  resistivity  for  the  flow  of 
hydrated  protons  (Q  cm),  t  the  thickness  of  membrane  (cm)  and 
Aactive  is  the  active  cell  area  (cm2).  Based  on  the  specific  fuel 
cell  plate  and  flowfield  design,  the  active  cell  area  is  only  a 
portion  of  the  stack  cross-sectional  area,  at  56%  [20].  The  other 
44%  of  the  area  are  used  to  accommodate  the  rods  that  hold  the 
stack,  the  manifolds  for  hydrogen  fuel,  oxidant  air  and  coolant, 
as  well  as  the  seals  surrounding  the  parameter  of  the  plate  and 
these  manifold  holes.  This  value  is  stack  design  dependent,  but 
representative  for  a  PEM  stack.  The  active  area  is  then  given  by 


Aactive  —  0.56  At 


(18) 


where  At  is  the  total  stack  cross-sectional  area.  While  the  other 
two  terms  in  Eq.  (17)  are  known  parameters  of  a  specific  cell, 
the  term  cm  is  difficult  to  describe  phenomenologically,  and  thus 
the  following  semi-empirical  expression  has  been  derived  [5] 

_  181. 6[1  +  0.03(0  +  0.062(7’ceii/303)2/2-5] 

™  "  [Vembrane  -  0.634  -  3i]  e3-25(rCe„-303/rcell)  (  9) 


where  again  i  is  the  current  density;  Amembrane  is  in  this  case  an 
adjustable  fitting  parameter  influenced  by  the  method  of  manu¬ 
facture  of  the  membrane,  and  a  function  of  the  relative  humidity 
and  stoichiometric  ratios  at  the  anode  and  cathode  and  of  the  age 
and  use  of  the  membrane.  The  parameter  is  usually  assigned  a 
value  between  10  and  20,  and  in  this  case  it  takes  the  value  of 
14.6  [5].  Eq.  (19)  is  developed  in  Ref.  [5]  in  a  semi-empirical 
manner,  based  on  mechanistic  equations  found  in  [22,23]  and 
experimental  data  at  RMC,  as  explained  in  [4] .  It  should  be  noted 
that  the  current  version  of  the  model  assumes  a  fully  humidified 
membrane;  a  membrane  hydration  module  is  currently  under 
development. 

2.1.4.  Concentration  overvoltage 

The  expression  used  in  the  model  for  the  concentration  over¬ 
voltage  is  given  by  [19] 

_  ^7"ceii  ^  f  A  1  \ 

Oconc  —  2^7  ln  (  1  j  J  (20) 

where  R  is  the  universal  gas  constant  (which  has  a  value  of 
8.3145  J  mol-1  K-1),  rcen  the  cell  temperature  (K),  F  Faraday’s 
constant  (which  has  a  value  of  96485.34  C  mol-1  electron-1), 
I  the  total  current  in  the  fuel  cell  (A)  and  /l  is  the  current  (A) 
at  which  the  hydrogen  fuel  is  used  up  at  a  rate  that  is  equal  to 
its  supply  rate,  hence  the  current  I  can  never  be  greater  than  /l. 
This  overvoltage  term  dominates  at  higher  current  densities. 

2.2.  Fuel  cell  gross  output 

The  gross  output  of  the  fuel  cell  stack  (W)  is  given  by 

fkgross, stack  =  7Ece]]/7ce]]  (21) 

where  nCQ\\  is  the  number  of  cells  in  the  stack.  The  number  of 
cells  is  27  for  this  study,  as  this  is  the  actual  number  of  cells  in 
the  physical  stack  on  which  testing  has  taken  place.  Note  that 
the  inclusion  of  the  cell  number  in  the  optimization  would  not 
yield  useful  results  since  there  is  no  penalty  for  increasing  the 
number  of  cells;  more  cells  would  simply  increase  the  output 
power  and  the  optimization  would  determine  the  upper  bound 
to  be  the  “optimal”  choice.  If  volumetric  and  system  weight 
considerations  were  taken  into  account  and  the  output  power 
decreased  with,  for  example,  increased  weight,  the  number  of 
cells  in  the  stack  could  be  included  in  the  optimization. 

2.2.1.  Heat  production 

Heat  is  generated  by  the  operation  of  the  fuel  cell  since  the 
enthalpy  that  is  not  converted  to  electrical  energy  will  instead  be 
converted  to  thermal  energy.  This  will  raise  the  temperature  of 
the  fuel  cell  beyond  its  operating  temperature  range,  and  must 
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be  addressed  using  a  cooling  system.  The  model  in  this  work 
assumes  that  the  passage  of  reactant  air  through  the  cell  will 
account  for  some  of  the  heat  transfer,  while  the  cooling  sys¬ 
tem  will  account  for  the  removal  of  the  remaining  excess  heat 
to  maintain  the  proper  cell  temperature.  A  fuel  cell  stack  will 
generate  the  following  amount  of  heat  during  operation: 

G total  —  (^max  —  ^cell)^cell  (22) 

where  Em ax  is  the  maximum  voltage  obtained  if  the  hydrogen 
‘caloric  value’,  i.e.  the  heating  value  or  enthalpy  of  formation 
were  transformed  into  electrical  energy,  and  is  given  by 

AH 

Emsix  —  —  =  1.48  V  (23) 

nf 

where  AH  is  the  change  in  enthalpy. 

Heat  will  be  lost  to  the  surroundings  and  to  the  air  stream 
through  three  heat  transfer  processes:  radiation,  convection  and 
evaporation.  Thus,  the  amount  of  heat  lost  to  the  air  stream  and 
surroundings  can  be  calculated  using  the  following  expression: 

G air, total  =  Grad  T  Gconv, total  T  Gevap  (24) 

The  amount  of  heat  emitted  through  radiation  by  the  fuel  cell 
stack  to  its  surroundings  is  given  by 

(2 rad  =  aAceu(rstack  -  T$)  (25) 

where  a  is  the  Stefan-Boltzmann  constant  (which  has  a  value 
of  5.67  xlO  8Wm  2  K  4),  rstack  is  the  average  temperature 
the  stack  (assumed  to  be  the  same  as  rceii),  Tq  is  the  ambient 
air  temperature,  and  Acen  is  the  area  of  the  radiative  stack  body, 
given  by 

Acell  =  0-85  X  2 A p;nds  +  0.70  X  1.3 Asides  (26) 


where  hc onv,top  and  hc onv, bottom  are  the  heat  transfer  coefficients 
due  to  convection  of  the  top  and  bottom  of  the  fuel  cell  channel, 
respectively.  The  total  heat  lost  to  convection  is  then  given  by 
the  sum  of  Eqs.  (27)  and  (29): 

Gconv, total  —  Gconv, walls  T  G conv,top+bottom  (20) 

The  model  calculates  the  amount  of  heat  lost  to  evapora¬ 
tion,  using  an  assumption  that  90%  of  the  water  produced  is 
evaporated.  Further  research  is  to  be  carried  out  to  confirm  this 
assumption.  The  amount  of  heat  required  to  convert  a  mass  of 
liquid  water  to  vapour  water  is  equal  to  the  product  of  the  latent 
heat  of  vaporization  /vap  and  the  mass  of  the  water.  An  expression 
for  the  amount  of  heat  lost  to  evaporation  is  given  as  [21] 

Gevap  =  ^vap  0.90mn2O, produced  (21) 

where  the  equation  for  water  production  was  given  in  Eq.  (5). 

2.3.  Parasitic  power  considerations 

Several  auxiliary  systems,  known  as  the  balance  of  plant 
(BOP),  are  necessary  for  the  correct  operation  of  a  fuel  cell 
system.  The  most  important  auxiliary  systems  are  the  air  com¬ 
pressor  necessary  to  pressurize  the  oxidant  air  to  the  correct  fuel 
cell  operating  pressure;  the  humidifier  to  guarantee  that  the  fuel 
cells  are  properly  humidified  for  optimal  performance;  and  the 
cooling  system  to  maintain  the  cell  temperature.  All  of  these 
auxiliary  systems  are  listed  and  a  representative  model  for  each 
is  described  in  the  following  sections. 

The  BOP  components  will  draw  power  from  that  produced 
by  the  fuel  cell,  thereby  reducing  the  overall  power  output.  The 
total  parasitic  power  of  the  fuel  cell  stack  is  modeled  by  the 
equation: 


where  the  decimal  terms  determine  the  actual  surface  area  of  the 
fuel  cell  stack.  In  this  model,  85%  of  the  two  end  plates  is  the 
effective  heat  radiating  area,  and  70%  of  the  surrounding  area 
of  the  stack  is  exposed.  The  geometry  of  the  stack  effectively 
increased  the  heat  radiating  area  by  30%,  leading  to  a  coefficient 
of  1 .3.  The  heat  loss  due  to  convection  will  occur  differently  for 
the  various  parts  of  the  fuel  cell  interior.  The  convection  lost  by 
the  vertical  walls  of  the  cell  is  given  by  [21] 

Gconv, walls  =  1  •  12  Asides/zconv(Tstack  —  ^b)  (27) 


where  hc onv  is  the  heat  transfer  coefficient  of  the  cell  channel 
sides  (W  m-2  K-4),  and  is  calculated  using  Eq.  (28)  [20].  The 
coefficient  1.15  also  incorporates  the  effective  side  area  of  the 
stack  based  on  the  non-smooth  stack  geometry.  The  heat  transfer 
coefficient  is  given  by 


h 


conv 


k  air  Nil 
l 


(28) 


where  ka ir  is  the  thermal  conductivity  of  air  (W  m  2  K  “4),  Nu 
the  dimensionless  Nusselt  number  and  /  is  the  stack  length.  The 
heat  transfer  due  to  convection  on  the  top  and  bottom  of  the 
channel  is  given  by  [21] 


Gconv,top+bottom  —  (h  conv, top  T  h  conv, bottom X^stack  ^i))  A  ends 


(29) 


Approximately,  80%  of  the  total  heat  produced  will  be  carried 
away  by  the  cooling  system;  the  rest  of  the  heat  will  be  removed 
by  the  reactant  air.  The  power  (W)  required  to  pump  the  coolant 
water  is  given  as  [20] 


Wcool, 


pump  — 


mn2o 

Ph2o 


^drop, coolant 


5F; 


actor 


Opump  7  motor 


(23) 


where  ^Factor  is  a  factor  of  safety  (assigned  a  value  of  1.5)  to 
account  for  any  pressure  losses  that  are  not  considered  explic¬ 
itly,  Pdrop, coolant  is  the  pressure  drop  of  the  coolant  through  the 
coolant  loop,  and  r]pump  and  ^motor  are  the  efficiencies  of  the 
pump  and  motor,  respectively.  The  assumptions  have  been  made 
in  the  model  that  there  is  a  pressure  drop  of  the  coolant  water 
through  the  system  of  1  x  104  Pa  and  that  the  efficiencies  for 
the  pump  and  pump  motor  are  60%  and  80%,  respectively.  The 
mass  of  water  (kg  s-1)  required  to  perform  the  heat  transfer  to 
the  water  coolant  is  [20] 

Gtotal  —  Gair, total  x 

WCOol,H20  =  - Tifi - ^77 - 7  (24) 

cPn2  o 'A  stack  ^  cool, H20, ini 


where  Cp^0  is  the  specific  heat  of  water  at  constant  pressure  and 
7cooi,H2o,in  is  the  temperature  of  the  coolant  water  at  the  inlet  of 
the  stack. 
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Operation  of  the  humidifier  pump  also  contributes  to  the 
power  loss.  The  total  humidifier  power  consumed  is  given  by 


The  system  energy  efficiency,  also  known  as  first  law  effi¬ 
ciency,  is  calculated  using  the  equation  [22] : 


w  _  7ft  H20, produced  D  ^Factor 

W humid, pump  —  ^clrop, humid 


Pwater 


0  pump  0  mot  or 


(35) 


The  model  implementations  then  assume  that  the  pressure  drop 
in  the  humidifier  to  be  1  x  105  Pa  and  that  the  pump  and  motor 
efficiencies  are  the  same  as  for  the  cooling  system.  The  amount 
of  water  necessary  to  fully  humidify  the  hydrogen  is  [24] 


8-937^H20  . 

m  H20,inH2  =  p -  Sat  OTH2  (36) 

''cell  P  h20 

where  rhu2  is  the  mass  flow  of  hydrogen  gas,  given  by 


^h2  =  Nu2MU2  (37) 

and  Mh2  is  the  molar  mass  of  hydrogen  (2.016  x  10-3 
kg  mol-1). 

Similarly,  the  amount  of  water  injected  into  the  air  stream  is 
given  by 

0-62 19/^0  . 

7ftH20,inair  —  „  sat  ft^air  (3o) 

“cell  —  Pu20 

where  mati  is  the  mass  flow  of  the  reactant  air,  given  by 


Tftair  —  N o2  37air  (39) 

where  Ma ir  is  the  molar  mass  of  air  (28.97  x  10-3  kg  mol-1). 
The  amount  of  water  passing  through  the  heat  exchanger  is  thus 


7ftH20,heat  exchanger  —  ^H20, produced  ^H20,in  H2  ^H20,in  air 

(40) 

Proof  of  the  connections  to  phenomenologically  based  theory 
for  Eqs.  (35)-(38)  has  not  yet  been  fully  established,  and  the 
equations  should  be  studied  in  detail  in  the  future  to  determine 
their  accuracies. 

A  compressor  is  required  only  for  the  cathode  side  of  the 
fuel  cell  since  the  hydrogen  gas  is  assumed  to  be  stored  in  a 
pressurized  container.  The  expression  for  the  compressor  power 
is  given  as  [19] 

Compressor  =  CpAm,J^  (  ( ^  ^  -  l)  (41) 

7? comp  \  V  ^air,atm  J  J 

where  c/?ajr  is  the  specific  heat  of  air  that  varies  with  temperature, 
y  is  the  isentropic  compression  ratio,  or  Cp/Cv,  and  is  calculated 
from  the  ratio  of  operation  pressure  PCeii  1°  atmospheric  pres¬ 
sure,  the  compressor  efficiency  is  assigned  a  value  of  60%,  and 
Pair, atm  is  atmospheric  pressure. 


2.4.  System  net  output  and  efficiency 

The  net  system  power  output  is  the  gross  output  power  of  the 
stack,  given  in  Eq.  (21),  subtracted  by  the  total  parasitic  power, 
calculated  in  Eq.  (32): 

fknet, system  —  fkgross, stack  —  ^net,  parasitic  (42) 


_  fknet,  system  rA^\ 

^system, energy  —  ^  ,  7  7 

Z^  in 7  v  react,  z  ft  react,  z  out 1 v  prod, ft  prod,  j 

where  Areact,z  is  the  molar  flow  (mol  s  )  and  /ireact,i  the  enthalpy 
(kJmol-1)  of  a  reactant  of  the  electrochemical  reaction  taking 
place  in  the  fuel  cell  stack  and  Nvmoj  is  the  molar  flow  (mol  s  ) 
and  /zprod,/  is  the  enthalpy  (kJ  mol  1 )  of  a  product.  The  inlet  molar 
mass  flow  of  H2  was  given  in  Eq.  (3).  The  enthalpy  terms  of  Eq. 
(43)  are  calculated  using  thermodynamic  tables  and  depend  only 
on  temperature  for  the  gases  (assumed  ideal)  and  on  the  temper¬ 
ature  and  enthalpy  of  formation  for  the  liquid  water  formed  at 
the  cathode. 

The  stoichiometry  of  hydrogen  in  this  study  is  set  to  be  1 . 1  and 
the  outlet  molar  flow  of  H2  (mol  s-1)  is  given  by  the  subtraction 
of  the  usage  from  the  supply: 


Afit2,out  —  ^H2,in 


7ft  cell 
IF 


1)  =  0.1 


7ft  cell 

IF 


(44) 


Again,  a  similar  expression  can  be  found  for  the  oxygen  at  the 
outlet  of  the  stack. 

Another  more  useful  metric  for  measuring  the  stack  perfor¬ 
mance  is  the  exergetic  efficiency  since  it  provides  a  true  analysis 
of  the  irreversibilities  of  the  system  that  the  energy  efficiency 
does  not  capture.  The  exergetic  efficiency  is  also  known  as  the 
second  law  efficiency.  The  equation  for  the  system  exergetic 
efficiency  is  given  by  [25] 


7? system, exergetic  — 


net, system 


y^in^react  .i  §react,z  X^out^prod  j'^prodj 

—  Gexit(l  —  To/  Tstack) 


(45) 


where  §react,z  is  the  exergy  (also  known  as  the  availability)  of 
a  reactant  i  flowing  into  the  system,  §prod j  is  the  exergy  of  a 
product  j  flowing  out  of  the  system,  Gexit  is  the  heat  leaving 
the  stack  (and  used  to  heat  the  air  and  release  the  hydrogen  if  a 
metal  hydride  hydrogen  storage  device,  instead  of  a  compressed 
hydrogen  tank,  is  used  as  the  onboard  fuel  storage  means),  Tfi 
is  the  ambient  air  temperature,  and  rstack  is  the  stack  operating 
temperature  (equal  to  rcen).  The  only  reactant  containing  exergy 
is  actually  the  hydrogen  stream;  the  product  hydrogen,  oxygen, 
nitrogen  and  water  exiting  the  system  will  all  possess  exergy. 
The  exergy  of  a  substance  (kJkmol-1)  is  given  by  [26] 

§  =  §KE  +  §PE  +  §TM  +  §CH  (46) 

where  §ke,  §pe,  §tm,  and  §ch  are  the  exergies  due  to  an  imbal¬ 
ance  of  the  substance  with  the  reference  kinetic  energy,  poten¬ 
tial  energy,  thermomechanical,  and  chemical  equilibrium  point, 
respectively.  The  first  two  terms  are  assumed  to  be  negligible, 
and  the  latter  two  terms  are  defined  by  [23] 

§tm  =  (h  -  ho)  -  To(s  -  So), 

§ch  =  7 (A6 jo  -  l1  joo)  (47) 

j 
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where  h ,  s,  and  /xy  denote  the  specific  enthalpy,  entropy,  and 
chemical  potential  of  substance  j ,  respectively,  at  a  given 
state,  and  xj  is  the  molar  fraction  of  substance  j.  The  sub¬ 
script  0  in  the  first  equation  indicates  that  the  property  is 
associated  with  the  reference  environmental  (restricted)  state 
where  pressure  and  temperature  are  at  equilibrium.  The  dou¬ 
ble  subscript  in  the  second  equation  refers  to  the  dead  state 
where  pressure,  temperature  and  chemical  potential  are  at 
equilibrium. 

3.  PEMFC  stack  and  empirical  parameters 

The  stack  used  for  the  study  is  based  on  the  Palcan  PC6- 
1200,  which  is  a  PEMFC  stack,  comprised  of  25  adjacent  cells 
in  a  conventional  plate-and-frame  arrangement  with  volumet¬ 
ric  dimensions  of  approximately  15  cm  x  19  cm  x  21  cm,  and 
its  total  active  area  is  96  cm2.  The  weight  of  the  stack  is  approx¬ 
imately  1.5  kg.  Inlet  and  exhaust  ports  for  fuel,  oxidant  and 
water  are  accessible  from  both  ends  of  the  stack  for  installation 
versatility  and  convenience.  The  designed  minimum  operating 
conditions  of  the  PC6  stack  are  pressures  slightly  above  atmo¬ 
spheric  (3psi  gauge,  or  0.21  bar  gauge)  and  a  temperature  of 
45  °C  average.  The  maximum  gross  power  output  of  the  stack 
(and  system)  is  1200  W.  The  empirical  coefficients  for  the  Mark 
V  PEMFC  stack  model,  given  in  the  last  row  of  the  three  entries 
in  Table  1 ,  are  used  in  this  study.  Extensive  stack  tests  on  the  Pal¬ 
can  PEMFC  stack  have  been  carried  out,  and  the  new  parameter 
models  resulting  from  these  tests  will  be  presented  in  subsequent 
work. 

4.  System  optimization 

The  design  of  a  PEMFC  system  often  involves  several  con¬ 
flicting  objectives.  In  this  case,  the  power  consumed  by  the 
auxiliary  systems  creates  several  conflicts  of  interest  when  the 
net  output  power  is  to  be  maximized.  The  operation  of  the  fuel 
cell  is  ameliorated  significantly  at  a  high  temperature  and  pres¬ 
sure;  however,  these  increases  would  result  in  higher  power 
losses  due  to  the  auxiliary  systems,  in  particular  due  to  the  cool¬ 
ing  system  and  the  air  compressor,  respectively.  It  is  because  of 
these  conflicting  interests  that  the  design  of  a  fuel  cell  system  is 
a  good  candidate  for  optimization. 

To  solve  the  optimization  problem  posed  by  the  conflict¬ 
ing  goals  in  the  fuel  cell  system  three  popular  and  successful 
optimization  techniques  have  been  chosen:  ( 1 )  simulated  anneal¬ 
ing  (SA);  (2)  genetic  algorithm  (GA);  (3)  sequential  quadratic 
programming  (SQP).  The  first  two  are  global  optimization  algo¬ 
rithms,  while  the  third  is  a  local  optimization  algorithm.  The 
global  optimization  algorithms  are  used  to  ensure  that  the  global 
minimum  is  determined  in  the  optimization.  The  local  optimiza¬ 
tion  algorithm  is  chosen  to  gain  some  insight  into  the  objective 
function  for  this  type  of  problem.  If  the  SQP  algorithm  does  not 
produce  results  similar  to  those  of  the  global  optimization  meth¬ 
ods,  then  it  is  clear  that  the  objective  function  is  not  unimodal 
and  smooth,  but  in  fact  contains  local  minima,  which  “trapped” 
the  local  optimization  algorithm. 


Design  variables 
(p,  T,  airStoich ) 


airStoich) 


Fig.  2.  Flow  diagram  of  the  optimization  program. 


To  solve  the  optimization  problem  the  fuel  cell  system  model 
described  in  the  previous  section  is  coupled  to  the  optimiza¬ 
tion  algorithms  described  above  as  shown  in  Fig.  2  for  the 
case  of  average  net  power  optimization.  The  coupling  between 
the  optimization  algorithm  and  the  model  is  achieved  through 
the  design  variables  and  the  objective  function  and  constraints. 
First,  the  optimization  algorithm  or  the  user  selects  the  initial 
value  for  the  design  variables  and  the  objective  function.  In  the 
case  of  global  methods,  the  initial  value  for  the  design  vari¬ 
ables  is  selected  randomly.  In  the  case  of  a  local  method,  the 
user  gives  the  initial  value  for  the  design  variables.  The  value 
of  the  design  variables  is  then  given  to  the  fuel  cell  system 
model.  The  fuel  cell  system  model  then  computes  the  per¬ 
formance  of  the  system,  and  from  those  results  computes  the 
value  of  the  objective  function  subject  to  the  design  constraints. 
With  this  information,  the  optimization  algorithm  chooses  a 
new  set  of  design  variables  that  can  potentially  increase  the 
system  performance.  This  process  is  repeated  until  a  conver¬ 
gence  criterion  is  satisfied  or  the  maximum  number  of  iter¬ 
ations  is  achieved,  depending  on  the  optimization  algorithm 
used. 

Finally,  if  a  local  optimization  algorithm  is  used  for  opti¬ 
mization,  the  gradients  are  computed  using  adaptive  forward 
differences.  In  this  case,  the  optimization  algorithm  automati¬ 
cally  calls  the  fuel  cell  system  model  with  a  small  perturbation 
in  each  one  of  the  design  variables  and  computes  the  numerical 
gradient. 

5.  Discussion  and  results 

The  fuel  cell  system  model  and  the  optimization  methods 
described  above  are  used  to  optimize  the  operating  conditions 
of  a  fuel  cell  system  for  two  different  applications:  vehicular  and 
stationary.  In  both  optimization  problems,  two  optimizations  are 
performed  with  different  objectives:  the  maximization  of  the  net 
power  and  of  the  system  exergetic  efficiency. 
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5.1.  Vehicular  application — a  multi- objective  optimization 
problem 

One  of  today’s  most  important  applications  of  PEM  fuel  cell 
system  is  to  serve  as  the  powerplant  for  an  electric  or  hybrid  vehi¬ 
cle.  In  this  application,  the  fuel  cell  system  needs  to  operate  over 
a  wide  range  of  power  requirements,  and  the  power  demands  on 
the  fuel  cell  system  change  with  varying  levels  of  acceleration 
and  velocity.  Therefore,  it  is  desired  to  obtain  the  fuel  cell  sys¬ 
tem  operating  conditions  that  yield  the  maximum  net  power  or 
exergetic  efficiency  over  varying  power  demands  to  maximize 
acceleration  and  velocity  or  vehicle  range  performance,  respec¬ 
tively.  The  fuel  cell  system  design  optimization  problem  is  then 
formulated  as  a  multi-objective  optimization  problem: 

N 

maximize  :  /(x,  i)  =  WjWj(x ,  //) 

w.r.t.x  (48) 

subject  to  :  xmax  >  x  >  xmin 
or 

N 

maximize  :  fix,  /)  =  uundx,  ij) 

w.r.t.x  ^  J  J  J  (49) 

subject  to  :  xmax  >  x  >  xmin 

where  W/(x,  ij)  is  the  net  power  and  ^(x,  ij)  the  system  exer¬ 
getic  efficiency  obtained  from  the  fuel  cell  system  at  a  current 
density  ij  for  objective  j ,  wj  the  weighting  of  objective  j,  and 
x  is  the  vector  of  design  variables.  The  net  power  output  and 
exergetic  efficiency  of  the  fuel  cell  system  is  obtained  using 
the  computational  model  described  above.  The  design  variables 
include  the  operating  temperature,  the  air  stoichiometry,  and 
the  operating  pressure  of  the  fuel  cell  stack.  The  upper  and 
lower  limits  of  these  design  variables  are  xmax=[373,  5,  15]  and 
Xmin  =  [338,  1.5,  1],  respectively.  If  only  two  objectives  were 
present,  it  would  be  possible  to  create  a  Pareto  curve  to  see  the 
trade-off  between  objectives;  in  this  case,  however,  the  objec¬ 
tives  are  the  maximization  of  the  average  net  power  and  average 
exergetic  efficiency  for  every  value  of  current  density  from  0.1 
to  1.3  A  cm-2.  The  multi-objective  optimization  problem  with 
different  weights  will  yield  different  optimal  operating  condi¬ 
tions.  The  weights  are  used  in  order  to  give  equal  importance  to 
the  power  or  efficiency  at  all  current  densities.  This  is  achieved 
by  using  the  weights  to  normalize  the  power  or  exergy  at  each 
current  density. 

The  maximization  of  the  net  power  is  solved  using  the  three- 
optimization  algorithms  described  above  and  the  results  are  pro¬ 
vided  in  Table  2.  The  solutions  of  the  three  methods  are  almost 


identical.  Only  the  solution  achieved  using  the  GA  has  a  small 
error  and  increasing  the  number  of  generations  in  the  algorithm 
will  reduce  this  error.  The  S  A  algorithm  stopped  when  a  conver¬ 
gence  criterion  of  three  successive  iterations  with  a  difference 
in  the  objective  function  smaller  than  1  x  10-6  was  satisfied. 
The  GA  algorithm  stopped  after  running  200  generations;  the 
convergence  criterion  was  not  specified.  Finally,  the  SQP  algo¬ 
rithm  stopped  when  the  module  of  the  gradient  of  the  objective 
function  reached  a  value  less  than  1  x  10~8.  In  principle,  the 
solution  of  the  SQP  is  dependent  on  the  given  initial  point.  The 
results  were  obtained  from  a  variety  of  different  initial  points 
and  the  same  solution  was  achieved  in  the  calculations. 

Due  to  the  different  nature  of  the  methods  and  the  different 
convergence  criteria  it  is  difficult  to  do  a  thorough  comparison 
of  the  computational  expense  of  these  methods.  Looking  at  the 
number  of  function  evaluations  and  the  required  computational 
time  in  Table  2  for  the  three  methods,  we  can  observe  that  the 
local  method  (SQP)  converges  to  the  solution  very  quickly,  and 
similar  results  were  obtained  when  the  initial  starting  point  was 
varied.  From  this  assessment,  it  can  be  observed  that  the  shape 
of  the  objective  function  is  likely  to  be  smooth  and  convex  (uni- 
modal)  since  the  local  method  converges  to  the  same  solution 
regardless  of  the  choice  for  the  initial  point. 

In  order  to  be  able  to  show  the  improvements  given  by  the  new 
operating  conditions,  the  performance  of  the  fuel  cell  system  is 
compared  to  the  performance  of  the  same  fuel  cell  at  nominal 
operating  conditions.  The  fuel  cell  system  running  at  the  nominal 
operating  condition  is  known  as  the  base  case  in  the  rest  of  the 
paper.  The  base  operating  conditions  include  a  temperature  of 
355  K,  hydrogen  and  air  stoichiometry  ratios  of  1.1  and  2.5, 
respectively,  and  a  pressure  of  3  bar.  Then,  performances  of  the 
fuel  cell  system  operating  at  the  calculated  optimal  operating 
conditions  and  at  base  operating  conditions  are  compared.  The 
differences  between  the  base  and  optimized  operating  conditions 
were:  3.2%  (temperature);  —40%  (air  stoichiometry);  —31.9% 
(pressure). 

Fig.  3  depicts  the  net  output  power  of  the  fuel  cell  system  at 
the  optimal  and  base  operating  conditions.  A  large  increase  in 
the  net  output  power  can  be  observed,  with  the  optimal  peak  at 
932.74  W  and  the  base  peak  at  728.36  W,  a  significant  increase 
of  21.9%.  This  increase  in  net  power  is  due  to  a  decrease  in  the 
necessary  power  of  the  auxiliary  devices  since  the  gross  output 
power  between  the  base  and  optimized  operating  conditions  is 
similar,  as  shown  in  Fig.  4.  The  optimization  has  also  moved 
the  point  at  which  maximum  net  power  is  achieved  to  a  higher 
current  density,  from  0.89  to  1 .007  A  cm-2,  a  change  of  1 1 .6%. 

At  the  base  operating  conditions  and  at  an  arbitrary  current 
density  of  0.76  A  cm-2,  the  consumed  powers  in  the  cooling 


Table  2 


Solution  of  the  multi-objective  problem  (net  power  optimization) 


Optimization  algorithm 

Objective  function: 
average  net  power  (W) 

Solution:  {T(K),  AirStoich,  P  (bar)} 

Number  of  function 

evaluations 

CPU  time  (s) 

Simulated  annealing 

663.68 

{366.84,  1.500,  2.043} 

3001 

55.1 

Genetic  algorithms 

663.80 

(366.41,  1.500,  2.011} 

4000 

615.8 

SQP 

663.68 

{366.84,  1.500,  2.043} 

100 

3.0 
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Current  Density  [A/cm2] 

Fig.  3 .  Net  output  power  of  the  fuel  cell  system  vs.  current  density  under  optimal 
and  base  operating  conditions. 

pump,  humidifier  pump  and  air  compressor  are  5.0  W,  5.0  W 
(minimum  allowed  power)  and  350.1  W,  respectively.  At  the 
optimal  operating  conditions  and  identical  current  density,  the 
cooling  pump,  humidifier  pump  and  air  compressor  consume 
powers  of  5.0,  5.0  and  130.4  W,  respectively.  This  represents 
decreases  of  0%,  0%,  and  62.8%  in  the  power  consumption  of 
each  parasitic  device.  The  parasitic  powers  of  the  cooling  and 
humidifier  pumps  are  unchanged  since  5.0  W  is  defined  as  the 
minimum  power  draw  of  the  devices,  and  both  the  base  case  and 
optimal  solution  show  that  this  is  more  than  sufficient  power  to 
fulfill  the  cooling  and  humidification  requirements.  The  large 
increase  in  net  power  is  mainly  due  to  the  massive  reduction  in 
the  power  consumed  by  the  air  compressor,  caused  by  a  large 
reduction  in  the  air  stoichiometric  ratio  and  operating  pressure. 
The  compressor  power  is  reduced  dramatically  at  all  current 
densities,  as  shown  in  Fig.  5. 

The  efficiency  of  the  fuel  cell  system  is  also  dramatically 
improved  by  the  optimization.  Fig.  6  illustrates  the  energy  effi¬ 


Fig.  4.  Gross  power  output  of  the  fuel  cell  stack  vs.  current  density  under  base 
optimal  and  base  operating  conditions. 


Fig.  5.  Consumed  power  by  the  air  compressor  vs.  current  density  under  optimal 
and  base  operating  conditions. 

ciency  as  a  function  of  current  density  for  both  the  optimal  and 
base  cases,  while  Fig.  7  depicts  the  exergetic  efficiency  as  a  func¬ 
tion  of  current  density.  Both  figures  demonstrate  that  significant 
efficiency  improvement  has  been  achieved.  The  peak  energy  effi¬ 
ciency  is  0.475  and  the  peak  exergetic  efficiency  is  0.700  (both 
at  0.06  A  cm-2)  for  the  optimized  operating  conditions,  while 
the  peak  energy  efficiency  is  0.425  (at  0.05  A  cm-2)  and  peak 
exergetic  efficiency  is  0.648  (at  0.06  A  cm-2)  for  the  base  case, 
for  percent  improvements  of  10.5%  and  7.4%,  respectively. 

It  is  also  possible  and  interesting  to  see  the  difference  in  per¬ 
formance  when  the  design  objective  is  to  maximize  the  system 
efficiency.  The  optimization  results  when  the  objective  function 
is  the  average  exergetic  efficiency  are  included  in  Table  3. 

In  this  different  design  optimization  problem,  the  maximum 
average  net  system  efficiency  obtained  was  0.515.  The  peak 
value  of  0.702  (at  0.11  A  cm-2)  is  slightly  higher  than  in  the 
previous  average  net  power  optimization  by  0.3%;  more  impor¬ 
tantly,  however,  the  current  density  at  which  the  peak  exergetic 


Fig.  6.  System  energy  efficiency  vs.  current  density  under  optimal  and  base 
operating  conditions. 
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Table  3 


Solution  of  the  multi-objective  problem  (exerg etic  efficiency  optimization) 


Optimization  algorithm 

Objective  function: 
average  exergetic  efficiency 

Solution:  {T (K),  AirStoich,  P  (bar)} 

Number  of  function 

evaluations 

CPU  time  (s) 

Simulated  annealing 

0.515 

{360.72,5.000,  1.036} 

3301 

23.4 

Genetic  algorithms 

0.515 

(360.72,5.000, 1.036} 

4000 

306.9 

SQP 

0.515 

(360.72,5.000,  1.036} 

189 

3.7 

Fig.  7.  System  exergetic  efficiency  vs.  current  density  under  optimal  and  base 
operating  conditions. 

efficiency  is  achieved  is  at  a  higher,  more  useful  level.  The  exer¬ 
getic  efficiency  as  a  function  of  current  density  is  shown  below 
in  Fig.  8.  The  peak  net  system  power  is  920. 16  (at  1 .01  A  cm-2), 
which  is  1 .3%  lower  than  in  the  case  of  average  net  system  power 
optimization.  An  optimization  of  the  exergetic  efficiency  of  a 
system  is  useful  because  it  can  directly  contribute  to  lower  ini¬ 
tial  investment  and  operating  costs.  The  shortage  of  maximum 
average  power  can  be  addressed  by  using  a  hybrid  power  system 


Fig.  8.  System  exergetic  efficiency  vs.  current  density  under  optimal  and  base 
operating  conditions  (exergetic  efficiency  optimization). 


with  a  battery  or  ultracapacitor  to  boost  the  system  power  when 
required. 

5.2.  Stationary  power  generation  application — a 
single-objective  optimization  problem 

Another  important  application  of  PEM  fuel  cell  systems  is  to 
serve  as  the  power  plant  for  a  mobile  or  stationary  power  unit 
to  provide  stable  power  supply.  In  this  case,  the  fuel  cell  sys¬ 
tem  is  to  be  operated  at  either  the  maximum  net  power  or  point 
of  maximum  exergetic  efficiency.  For  the  simple  power  gener¬ 
ation  application,  it  is  possible  to  maintain  the  fuel  cell  system 
at  a  single  operating  point.  Therefore,  the  fuel  cell  operating 
conditions  and  the  current  density  that  yields  the  maximum  net 
power  or  maximum  exergetic  efficiency  are  to  be  obtained.  The 
optimization  problem  can  be  formulated  as 

maximize  :  fix)  =  fk(x) 

wrt- x  (50) 

subject  to  :  xmax  >  x  >  xmin 

maximize  :  /(x)  =  ^(x) 

w-r-tx  (51) 

subject  to  :  xmax  >  x  >  xmin 

where  W(x)  is  the  net  power;  rj(x)  is  the  system  exergetic  effi¬ 
ciency;  x  is  the  design  variable  vector  that  includes  the  operating 
temperature,  air  stoichiometry,  operating  pressure  of  the  fuel  cell 
stack,  and  current  density  at  which  the  fuel  cell  system  will  oper¬ 
ate  to  achieve  maximum  net  power.  The  upper  and  lower  limits 
of  the  design  variables  are  defined  as  xmin  =  [338,  1.5,  1,  0.1] 
and  xmax  =  [373,  5,  15,  1.3]. 

This  problem  is  also  solved  using  the  three- optimization  algo¬ 
rithms  described  and  the  solutions  are  presented  in  Table  4. 
Furthermore,  two  initial  design  points  were  used  to  solve  the 
problem  using  the  local  method.  For  SQP-1,  xo  =  [353,  2.5,  3, 
0.75]  and  for  SQP-2,  x0  =  [353,  2.5,  1,  0.75]. 

In  this  case,  there  are  two  different  solutions  with  widely 
varying  operating  conditions  found  by  the  global  and  local  algo¬ 
rithms,  although  there  is  again  a  small  discrepancy  observed 
between  the  solution  obtained  using  the  GA  and  SA  optimiza¬ 
tion  methods,  even  though  they  in  effect  found  the  same  solution. 
As  in  the  previous  case,  an  increase  in  the  number  of  generations 
or  the  population  in  GA  would  likely  result  in  a  better  result  at 
the  expense  of  a  longer  computation  time.  The  existence  of  two 
distinct  solutions  suggests  that  the  objective  function  is  non- 
unimodal.  One  solution,  located  by  the  local  SQP-1  algorithm, 
is  at  a  very  high  temperature  and  high  pressure  with  low  air  sto¬ 
ichiometry  (the  minimum  value  is  reached);  the  other  solution, 
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Table  4 

Optimal  solution  of  the  single-objective  problem  (net  power  optimization) 
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Optimization  algorithm 

Objective  function: 
maximum  net  power  (W) 

Solution:  {T  (K),  AirStoich,  P  (bar),  i  (A/cm2)} 

Number  of  function 

evaluations 

CPU  time  (s) 

Simulated  annealing 

935.12 

{353.21,5.000,  1.028,0.984} 

4401 

7.8 

Genetic  algorithms 

934.65 

(353.33,4.578,  1.029,0.982} 

4000 

92.6 

SQP-1 

933.53 

(368.17,  1.500,  1.952,  1.007} 

140 

0.2 

SQP-2 

935.12 

(353.21,5.000,  1.028,0.984} 

244 

0.3 

found  by  both  global  algorithms  and  the  SQP  algorithm  with  a 
different  initial  starting  point,  is  at  a  much  lower  temperature 
and  pressure,  but  with  a  high  air  stoichiometry  (the  maximum 
value  is  reached). 

The  operating  conditions  for  the  solution  located  by  the 
SQP  algorithm  (SQP-1)  are  similar  to  the  solution  of  the  multi¬ 
objective  function  problem,  although  a  distinction  between  the 
two  results  is  obvious.  The  optimal  temperature  is  greater  for  the 
stationary  application  than  for  the  vehicular  application,  while 
the  optimal  operating  pressure  is  less.  The  optimal  air  stoichiom¬ 
etry  is  the  same  for  both  applications,  at  the  lower  bound  of  the 
design  variable  limits.  The  difference  between  the  maximized 
net  power  curves  shows  that  a  slight  compromise  was  required 
between  the  performances  at  high  and  low  current  densities  for 
the  multi-objective  optimization  problem:  the  peak  net  system 
power  for  the  multi- objective  problem  is  0.08%  less  than  the 
peak  net  system  power  determined  for  the  single-objective  prob¬ 
lem. 

Fig.  9  shows  a  curve  of  the  maximum  net  output  power  with 
respect  to  the  current  density  at  the  base  and  optimal  operating 
conditions  for  the  solution  found  by  the  SQP-1  algorithm.  The 
peak  of  the  optimized  net  system  power  curve  is  denoted  by 
the  small  (black)  circle  at  the  maximum  value.  It  is  interesting 
to  note  that  the  maximum  net  power  is  achieved  at  the  same 
current  density,  1 .007  A  cm-2,  for  which  the  maximum  average 
net  power  was  determined  in  the  multi- objective  problem.  The 
gross  power  curves  are  shown  in  Fig.  10. 


Fig.  9.  Fuel  cell  system  net  output  power  vs.  current  density  with  peak  power 
shown  (SQP-1  algorithm). 


Fig.  10.  Gross  power  output  of  the  fuel  cell  stack  vs.  current  density  under  base 
optimal  and  base  operating  conditions  (SQP-1  algorithm). 

Figs.  1 1  and  12  depict  the  net  and  gross  system  power  for  the 
global  algorithms’  solution.  The  net  power  curve  is  essentially 
the  same  as  that  for  the  SQP- 1  algorithm,  with  slight  differences 
in  peak  net  power  values  and  current  density  at  which  the  peak 
net  power  is  achieved.  However,  the  gross  power  curve  is  quite 
different  showing  that  the  gross  power  of  the  fuel  cell  stack 
has  been  reduced  significantly,  due  to  the  large  reduction  on 
the  operating  pressure  and  temperature.  Even  though  the  gross 


Fig.  11.  Fuel  cell  system  net  output  power  vs.  current  density  with  peak  power 
shown  (global  algorithm). 
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Table  5 

Optimal  solution  of  the  single-objective  problem  (system  exergetic  efficiency) 


Optimization  algorithm 

Objective  function:  maximum 
exergetic  efficiency  (W) 

Solution:  {T  (K),  AirStoich,  P  (bar),  i  (A/cm2)} 

Number  of  function 

evaluations 

CPU  time  (s) 

Simulated  annealing 

0.729 

{363.13,5.000, 1.158,0.100} 

4401 

5.3 

Genetic  algorithms 

0.729 

{363.31,4.997,  1.159,  0.100} 

4000 

64.8 

SQP 

0.722 

{361.23,4.623, 1.175,0.100} 

506 

0.8 

Fig.  12.  Gross  power  output  of  the  fuel  cell  stack  vs.  current  density  under  base  Fig.  14.  System  exergetic  efficiency  vs.  current  density  under  optimal  and  base 
optimal  and  base  operating  conditions  (global  algorithm).  operating  conditions  for  exergetic  efficiency  optimization. 


power  is  reduced,  the  increase  in  net  power  can  be  explained 
by  looking  at  Fig.  13,  which  shows  the  compressor  power  for 
the  global  algorithms’  solution.  When  contrasted  with  the  com¬ 
pressor  power  curve  of  Fig.  5,  a  dramatic  reduction  in  power 
can  be  easily  observed.  This  reduction  in  compressor  power  is 
mainly  due  to  the  reduction  of  the  pressure  ratio  in  Eq.  (41).  As 
the  pressure  ratio  tends  to  one  the  compressor  power  goes  to 
zero.  This  proves  that  there  is  a  large  penalty  for  increasing  the 
cathode  pressure  of  the  fuel  cell  stack. 
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Fig.  13.  Consumed  power  by  the  air  compressor  vs.  current  density  under  opti¬ 
mal  and  base  operating  conditions  (global  algorithm). 


As  in  the  previous  multi-objective  optimization  problem,  the 
system  exergetic  efficiency  is  optimized.  The  results  of  the  opti¬ 
mization  process  are  depicted  below  in  Table  5. 

The  maximum  peak  exergetic  efficiency  was  0.729,  which 
is  an  increase  from  the  previous  single-objective  of  4.0%  and 
4.3%  for  the  low  and  high  stoichiometry  solutions,  respectively. 
The  peak  of  the  net  system  power  fell  below  that  of  the  previ¬ 
ously  optimized  average  to  868.8 1  W  (at  0.96  A  cm-2),  although 
it  remained  higher  than  that  of  the  base  operating  case.  The 
exergetic  efficiency  curve  as  a  function  of  the  current  density 
is  shown  in  Fig.  14.  As  stated  previously,  the  maximization  of 
the  system  exergetic  efficiency  rather  than  the  net  system  power 
can  reduce  the  overall  costs  of  the  system,  and  the  shortfall  in 
power  can  be  mitigated  through  the  introduction  of  a  battery  or 
ultracapacitor  in  a  hybridization  scheme. 

6.  Conclusions 

In  this  paper,  a  methodology  to  obtain  the  optimal  operating 
conditions  for  a  fuel  cell  system  is  outlined.  In  order  to  show 
the  validity  of  the  approach,  the  operating  conditions  for  a  fuel 
cell  system  employed  in  two  different  applications  is  shown. 
The  new  operating  conditions  improve  the  net  output  power  or 
the  exergetic  efficiency  with  respect  to  the  nominal  or  common 
operating  conditions  of  the  system  for  both  applications,  demon¬ 
strating  the  usefulness  of  the  approach. 

During  the  study,  global  and  local  optimization  algorithms 
were  compared.  For  the  multi-objective  problem,  similar  results 
are  obtained  using  both  methods.  This  allows  us  to  conclude  that 
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this  optimization  problem  has  an  objective  function  that  is  both 
convex  and  unimodal.  However,  the  addition  of  design  variables 
and  constraints  may  change  the  nature  of  the  problem,  and  the 
use  of  local  methods  may  not  be  possible.  In  this  particular  case, 
the  use  of  a  local  optimization  technique  is  preferred  since  its 
application  results  in  the  same  solution  as  the  global  methods 
using  less  computational  time.  For  the  single- objective  optimiza¬ 
tion  problem,  the  differing  results  of  the  two  methods  suggest 
that  the  objective  function  is  in  this  case  non-unimodal  with  at 
least  one  local  maxima  in  addition  to  the  global  maximum  value. 
Furthermore,  the  benefit  of  using  both  method  types  is  demon¬ 
strated  by  the  large  variance  in  operating  conditions  revealed  by 
the  solutions:  in  practice,  one  maximum,  although  at  a  slightly 
lower  value,  may  actually  be  the  superior  choice.  This  interesting 
result  suggests  that  it  is  prudent  to  continue  to  employ  both  types 
of  algorithms  in  tandem  during  any  optimization,  if  possible. 

This  work  reveals  a  large  number  of  possible  avenues  of 
future  research.  The  existing  fuel  cell  stack  and  system  mod¬ 
els  would  benefit  from  fewer  assumptions  and  simplifications; 
for  example,  the  assumed  constant  pressure  drop  in  the  humidi¬ 
fier  and  in  the  cooling  channels  and  fully  humidified  membrane. 
A  new  optimization  could  be  performed  that  takes  into  account 
not  only  the  operating  conditions  as  design  variables  but  also 
the  geometric  parameters  of  the  fuel  cell.  Indeed,  a  combina¬ 
tion  of  operating  parameters  with  physico-chemical  parameters 
as  design  variables  would  result  in  a  comprehensive  optimiza¬ 
tion  of  a  fuel  cell  system.  Also,  the  objective  function  can  be 
modified  to  reflect  the  operating  conditions  of  an  actual  fuel 
cell  application:  the  weighting  factors,  set  to  1  for  this  study, 
could  be  matched  to  the  particular  power  cycle  of  the  applica¬ 
tion,  increasing  the  accuracy  of  the  model.  Subsequent  papers 
will  reflect  this  ongoing  attempt  to  increase  the  sophistication 
and  accuracy  of  the  stack  and  system  models. 
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